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Abstract 

We consider Schrddinger operators on possibly noncompact Rie- 
mannian manifolds, acting on sections in vector bundles, with locally 
square integrable potentials whose negative part is in the underlying 
Kato class. Using path integral methods, we prove that under geodesic 
completeness these differential operators are essentially self-adjoint on 
C^°, and that the corresponding operator closures are semibounded 
from below. These results apply to nonrelativistic Pauli-Dirac opera- 
tors that describe the energy of Hydrogen type atoms on Riemannian 
3-manifolds. 



1 Introduction 

A classical result from B. Simon's seminal paper [27] states that a Schrddinger 
operator of the form —A + V in the Euclidean space M m , with V: W 71 — > M 
a locally square integrable potential, is essentially self-adjoint on C[f(lR m ), if 
the negative part of V is in the Kato class JC(W m ). Note here that this fact is 
closely related to quantum physics, in the sense that the Coulomb potential 
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V(x) = —l/\x\ is in the above class. Having in mind that all of the above 
data can be defined on any Riemannian manifold, we are interested in the 
following question in this paper: 

To what extent can Simon's result be extended to Schrddinger type 
operators acting on sections in vector bundles over possibly non- 
compact Riemannian manifolds? 

Apart from a pure academic interest, this question is also particularly mo- 
tivated by the observation that it is possible to model [16, 9] nonrelativistic 
atomic Hamiltonians on any nonparabolic Riemannian 3-manifold (which 
have to be spin c , if the electron's spin is taken into account; in particular, 
the vector-valued case becomes particularly interesting from this point of 
view, see Section 3 below). This abstraction is desirable from the physics 
point of view, since one would like to understand deeply which properties 
of the Euclidean space M 3 actually guarantee certain spectral properties of 
quantum systems, or other important results such as the stability of mat- 
ter [23]. In these situations, the corresponding potential terms are always 
locally square integrable, and with some control on the underlying Rieman- 
nian structure, their negative parts are in the underlying Kato class, so that 
we basically are in the initial situation. 

Before we can formulate our main result, we have to introduce some 
notation: 

Let M denote a smooth connected Riemannian manifold without bound- 
ary. The geodesic distance on M will be written as d(x,y), and K r (x) will 
stand for the open geodesic ball with radius r around x, and 

(0, oo) x M x M — > (0, oo), (t, x, y) i — > p(t, x, y) 

will stand for the minimal positive heat kernel on M. 

If F — > M is a smooth Hermitian vector bundle, then, abusing the no- 
tation in the usual way, ||«|| x stands for the norm and the operator norm 
corresponding to (•, m) x on each (finite-dimensional) fiber F x , and the scalar 
product and norm corresponding to the Hilbert space I\2(M, F) will be writ- 
ten as (•,•) and ||«||, respectively, that is, 

</l,/2>= / (h(x)J 2 (x)) x VO\(dx), ||/|| 2 = I ||/(x)||ivol(dx). (1) 
JM JM 

We will also write 

:M->- [0,oo), |*| (x) := \\*(x)\\ x 
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for any section $ in F or in End(F). If F — > M is a second bundle as above 
and if 

P:r C oo(M,F)^r C go(M,F) 

is a linear differential operator, then we denote with P' the formal adjoint 
of P with respect to (1). In particular, the Laplace-Beltrami operator on M 
is given in this sense as —A = d*d. The symbol V™ will denote the Levi- 
Civita connection, and if nothing else is said, the (co-)tangent bundle of M 
will be equipped with the Hermitian structure corresponding to the under- 
lying Riemannian metric of M. These data will be implicitely complexified, 
whenever necessary. 

Let E — > M be a smooth Hermitian vector bundle, let V be a Hermitian 
covariant derivative in E and let V: M — > End(P) be a potential, that is, V 
is a measurable section in End(P) such that V(x) : E x — > E x is self-adjoint 
for almost every (a.e.) x G M. Furthermore, let /C(M) denote the class of 
Kato functions 1 on M. Our main result reads as follows: 

Theorem 1.1. Let M be geodesically complete, let \ V\ G Lf oc (M) and assume 
that V admits a decomposition V — V± — V2 into potentials V 3 , > with 
I V2 1 G /C(M). Then the operator V^V/2 + V is essentially self-adjoint on 
Tcg° (M, E) and its closure is semibounded from below. 

Note that the decomposition V — V± — V2 into nonnegative potentials 
need not be the canonic one given by V = V + — V~, which can be defined 
through the fiberwise spectral calculus of E. 

Before we explain the strategy of the proof of Theorem 1.1, some remarks 
are in order: 

Remark 1.2. (a) Theorem 1.1 is disjoint from the various results on essen- 
tial self-adjointness for operators of the form V^V/2 + V that have been 
obtained in [1]. The point here is that, in general, Kato potentials need not 
satisfy the inequality (2.2) from [1], i.e., for every compact K C M there are 
numbers < Or- < 1, > such that 

(J \\V2(x)\\l \u(x)\\ol(dx)y /2 < a J r||A«||+&* r H ( 2 ) 

for any u G Cjf (M). As an example, which also illustrates the kind of 
singularties of V 2 we can deal with, consider V^ix) := |x|~ 2 | log(|a;|)| _5 
in the Euclidean space R m , with 5 > 0. If m > 3, then V 2 {5) G /C(M m ) 
if and only if 5 > 1 ([4, p.8]), but for example, one has V 2 (2) ^ Lf oc (R 3 ), 

1 see Section 3.5 for the definition of K.(M) and for criteria for functions to be in JC{M) 



3 



whereas (2) clearly implies local square integrability. However, it should be 
noted that the main strength of the results of [1] lies in the fact that the 
authors have considered arbitrary first order elliptic differential operators 
instead of V. It would certainly be an interesting problem to see to what 
extent our probabilistic techniques below can be extended to cover the latter 
situation, which has first been considered in [22]. 

(b) Of course, taking E = M x C and V = d + i(5 with (5 E ^(M), we can 
deal with smooth magnetic potentials within our framework. In this scalar 
situation, the analogue of Theorem 1.1 can be easily deduced from (a slight 
variation of) Theorem 1 in [13], where the authors can even allow magnetic 
potentials with possibly strong local singularities. We refer the reader to [17] 
for the scalar situation in Euclidean space. 

Let us now explain the strategy (which is partially motivated by [27] 
and [13]) of the proof of Theorem 1.1, which is given in full detail in the 
following Section 3.5. To this end, we assume for the rest of this section that 
V is as in Theorem 1.1. Then by the main result of [14], it is always possible to 
define the form sum Hy corresponding to the Friedrichs realization of V^V/2 
and V without any additional assumptions on M (see Theorem 2.6 below). 
The main advantage of this observation is that, unlike in usual essential- self- 
adjointness proofs, instead of directly proving that V*V/2 + V is essentially 
self- adjoint on Tcg°(M, E), we will prove that the latter space is an operator 
core for Hy (this is the content of Theorem 2.14; Theorem 1.1 itself follows 
directly from the latter result, which is summarized in Corollary 2.15). In 
particular, we will use the full spectral calculus given by Hy. 

Having said this, the first step in the proof of this operator core property 
will be to deduce the following smoothing property (see Proposition 2.11 
below): 



where B(x) is a Brownian motion starting in x with lifetime ((x), where 



For any t > one has e~ tHv r L2 (M, E) C r L g> c (M, E). 



(3) 



This result will be derived from the path integral formula 



e- tH -f(x) = E [l {t< « x)} y t * IT' f(B t (x 



■))], 



(4) 




is the corresponding stochastic parallel transport with respect to V, //\ 
//t'* its inverse, and where 



■x,— 1 
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% x : E, 
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is the path ordered exponential 2 



% x - 1 



/ //r 1 nBsA^)n---//vnB Sk (x))//ids 1 ...ds k 



(5) 



(details on these processes and on formula (4), which is one of the main 
results of [15], are included in the following section). Again, (3) and (4) are 
valid without any additional assumptions on M. 

Remark 1.3. Note that it is not possible to deduce (3) directly by Sobolev 
embedding theorems for dimM > 3, which is the main motivation for the 
introduction of path integral techniques in this context. 

In a next step, we will use finite speed propagation methods to deduce 
the following result: 



To be precise, we will actually prove a Davies-Gaffney inequality (see 
Proposition 2.13) for approximations of Hy and use the fact that this in- 
equality always implies (is in fact equivalent) to finite speed of propagation 
by the results of [3]. Then one can use a variant of Chernoff's theorem (see 
Lemma B.l) to deduce (6). The fact that we use finite speed propagation 
methods in this context has been particularly motivated by the scalar situa- 
tion that has been considered in [13], where the authors apply this method 
in a similar way. As has been noted in [13], this technique avoids the usage 
of second order cut-off functions, which do not seem to be available without 
additional control on the underlying Riemannian structure. 

As a next step one can combine (6) with (3) to deduce the following fact: 



The set D(H V ) H IY- (M, E) n |/| / has a compact support j (7) 
is an operator core for H V) if M is geodesically complete. 

Then, we shall use the self-adjointness of H v to deduce that the elements 
/ of the set (7) satisfy V f V/ G r L2 (M,£). Finally, if M is geodesically 

2 Here, tAk — {0 < Si < ■ ■ ■ < Sk < t} C M fc denotes the i-scaled fc-simplex for any 
k e N, t > 0. 




(6) 
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complete we can use a (local) result on Friedrichs mollifiers to prove that 
r C g=(M, E) is an operator core for H v , by showing that r C g°(M, E) is dense 
in (7) with respect to the graph norm corresponding to Hy. 

This paper is organized as follows: In Section 2, we first recall some facts 
about Kato potentials. The rest of Section 2 is completely devoted to the 
proof of Theorem 1.1. In Section 3, we apply Theorem 1.1 in the context of 
Hydrogen type problems on Riemannian 3-manifolds, which was originally 
the main motivation for this paper. Finally, in the appendix, we have in- 
cluded a fact about Friedrichs mollifiers, an abstract variant of Chernoff's 
finite speed of propagation theorem on vector bundles, and some facts about 
path ordered exponentials that we will need in our probabilistic considera- 
tions. 

2 Kato potentials and the proof of Theo- 
rem 1.1 

Let us first clarify that in this section, 

M will always be a smooth connected Riemannian manifold with- 
out boundary, E — >■ M a smooth Hermitian vector bundle, V a 
Hermitian covariant derivative in E, and V: M — > End(.E') a 
potential. 

By the usual abuse of notation, we will denote the quadratic form corre- 
sponding to a symmetric sequilinear form in some Hilbert space with the 
same symbol. The symbol Ho stands for the Friedrichs realization of V^V/2, 
that is, Hq is the nonnegative self-adjoint operator in T\_2(M,E) which cor- 
responds to the closure q Ho of the quadratic form given by the symmetric 
nonnegative operator V t V/2, defined initially on r c °°(M, E). Note the well- 
known: 

Remark 2.1. If M is geodesically complete, then one has 

D(<te )= {/|/ er L2 (M,£),v/ er L2 (M,£®T*M)}, (8) 
QH (f,h) = l I (Vf(x),Vh(x)) x vol(dx), 

1 J M 

and r C g=(M, E) is an operator core for Hq, and one has 

D(H )= {/|/,V t V/er L2 (M, J E?)}, ifo/^V+V/. (9) 
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Next, we remark that V defines a quadratic form in T L 2(M, E) by setting 

D(gy)= {f\feT L ,(M,E), (Vf, f) g L 1 (M)|, 



5V(/)=/ (V(x)/(x),/(x)) xV ol(dx). (10) 

We will often require a global Kato assumption on some negative part 
of V. Before recalling some facts on Kato functions, let us first introduce 
some notation: Let ^# := (Q, J^", P) be a filtered probability space which 
satisfies the usual assumptions. We assume that jjt is chosen in a way such 
that carries an appropriate family of Brownian motions 

B(x) : [0, C(a:)) x Q — > M, x G M, 

where ((x) : — )■ [0, oo] is the lifetime of We will freely use the fact 

¥{B t (x) £ N,t < ((x)} = / pt(x, y)vol(dy) for any measurable N C M 

Jn 

in the following. 

Now a measurable function w : M — > C is said to be in the Kato class 
K{M) of M, if 



lim sup E 



1 {*<C(*)} \w(B 3 (x))\ds 



0, which is equivalent to (11) 



lim sup / / p(s, x,y) \w(y)\vo\(dy)ds = 0. (12) 
*^°+ xeM Jo Jm 

The local Kato class K,\ oc {M) is defined in the obvious way, 

/Ci oc (M) := lx-iu G /C(M) for all compact AT C m} d /C(M), 

and generally, /Ci oc (M) may depend on the Riemannian structure of M. 
For future reference, we note: 

Lemma 2.2. (a) One has K,{M) C Lj^M) and L°°(M) C K(M). 

(b) For any w G l_i oc (M) and a.e. x G M one has 

¥{w(B.(x))eLl oc [0,ax))} = l. (13) 

(c) For any w G K,i oc (M) and all x G M one has 

¥{w(B.(x))eLl oc [0,ax))} = l. 
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(d) For any w G K(M), t > 0, one has 



sup E 

x&M 



{t<a*)} e 



J*\w(B 3 (x))\ds 



< 00. (14) 



Proof. Part (a) is an elementary result which is included in [14], and the 
parts (b)-(d) are included in Prop. 2.4 and Prop. 2.5 in [15]. ■ 

Let us now point out that with some control on the Riemannian struc- 
ture of M, one can easily produce a large class of Kato functions. To this 
end, we first note the following highly nontrivial self-improvement result of 
on-diagonal upper estimates for p(t,x,y), which will be very useful in the 
following: 

Theorem 2.3. Assume that there is a C > and a T e (0, oo] such that 

C 

sup p(t, x, x) < tdimM/2 for allO<t<T. 
Then there are C\, C2 > such that 

sup pit, x, y) < -Jh—e-^l^ for allO<t<T. 

The reader may find a proof of this result in [12] (see Theorem 1.1 therein 
for a more general result). Theorem 2.3 easily implies: 

Proposition 2.4. Assume that there is a C > and a T > such that for 
all < t < T one has 

C 

supp(t,x,a;)<^--^. (15) 

Then for any p such that p > 1 if dim M = 1, and p > dim M/2 if dim M > 
2, one has 

L P (M) + L°°(M) C /C(M). (16) 

Proof. Indeed, Theorem 2.3 implies the existence of a C > such that for 
all < t < T one has 3 

C 

sup p{t, x, y) < 



x,y£M t ' 

Now we can directly apply Proposition 2.8 in [14] (the corresponding proof 
is elementary and essentially only uses Holder's inequality). ■ 



3 Of course this inequality can also be deduced with an elementary argument. 
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Remark 2.5. The estimate (15) is satisfied, if M is geodesically complete 
with Ricci curvature bounded from below and a positive injectivity radius 
(see for example [21], p. 110). We refer the reader to [14] and the references 
therein for several other aspects of Kato potentials. 

The following result is also included in [14]. It shows that, remarkably, 
one can always define the form sum of H and V under the following very 
weak assumptions on V: 

Theorem 2.6. Let V be such that there is a decomposition V — V\ — Vi into 
potentials Vj > with \Vi\ G \-l oc (M) and \V 2 \ G /C(M). Then one has 

V(q Ho + qv) = DWnDfe), (17) 

and qHo+Qv ^ a densely defined, closed and semibounded from below quadratic 
form in T L 2 (M, E) . 

In the situation of Theorem 2.6, the form sum Hq + V will be denoted 
with H v , that is, H v is the self-adjoint semibounded from below operator 
corresponding to qu + qv- 

Remark 2.7. In the situation of Theorem 2.6, assume that M is geodesically 
complete. Then Proposition 2.14 in [14] states that T C g= (M, E) is a form core 
fortfy. 

Let us add the following simple observation: 

Lemma 2.8. Let \V\ € Lf oc (M) and assume that there is a decomposition 
V — Vi — V2 into potentials Vj > with | V2 1 G K{M). Furthermore, let 
Hv,min denote the operator V^V /2 + V with domain of definition Tc^(M, E) , 

and let Hy m \^ := Hy^mm- Then one has -ffymm C Hy. 

Proof. Since Hy is closed, it is sufficient to prove Hy tm[n C Hy. But if 
/ G F C ~(M,E), h G D(q Hv ), then / G D(q Hv ) and we have 

QH v (f,h) = ±(tfVf,h) + <yf,h), (18) 

so / G D(H V ) and H v f = |V f V/ + Vf. ■ 

As we have already remarked in the introduction, an essential step in the 
proof of Theorem 1.1 will be to deduce an L 2 -w Lj^ c smoothing property of 
the Schrodinger semigroup 

(e- % ) f >oCif(r L2 (M,£)), 
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which will be deduced from a path integral formula for e~ tHv . In order to 
formulate the latter formula in our geometric context, for any t > the 
stochastic parallel transport with respect to (B(x), V) will be written as a 
pathwise unitary map 

// x t : E x — > E Bt(x) , defined in {t < £(x)} C ft. 

Now Theorem 2.11 in [15] states the following Feynman-Kac type path 
integral formula: 

Theorem 2.9. In the situation of Theorem 2.6, for a.e. x G M, there is a 
unique process 

f x : [0,C(a:)) x ft — >• End(£ x ) 

which satisfies 

c\i /x / -, \ 

= -^( //r 1 in ), n = i (19) 

pathwise in the weak sense, and for any f G T L 2(M, £ > 0, a.e. x G M 
one /ias 

e-^/(:r) = E [l {t<c(se ) } ^ //r-" 1 f(B t (x))} . (20) 

Remark 2.10. The set of x for which i /x exists is, by definition, equal to 
the set x for which one has (13) for w = \V\, and if x is in this set, then the 
asserted formula (5) from the introduction follows from Lemma C.l. 

We will use (20) to deduce: 

Proposition 2.11. In the situation of Theorem 2.6, one has 

T L 2(M,E)] C r L c° c (M,E) for any t > 0. (21) 



e -tH v 



Remark 2.12. Note that Lemma 2.2, Theorem 2.6, Theorem 2.9, and Propo- 
sition 2.11 are all valid without any further assumptions on the Riemannian 
structure of M. 

Proof of Proposition 2.11. We define scalar potentials Vj : M — >■ [0, oo), v : 
M -> R by 

fi(«) := mincr(Vi(«)), v 2 {») := maxa(V 2 (»)), f(«) := fi(«) - v 2 (»)- 
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Let x be such that (13) holds for w — \V\\ and w = |V 2 |. Then V x exists, 
and V > v 1, Lemma C.l and — v < v 2 imply 

II^ILl {t <CW}<e- /ot,,(Bs(a:))dSl {*<C(,)} 

< eti V2{Bs{x))ds l {t<c(x)} P-a.s. for any t > 0, 

so that for any t > one has 

I|e [i { t<c(*»*r nr 1 f(B t (x))]\\ x 

^E^^e/o-^^ll/^^))!!^^ 




!{*<«.)} \\f(B t (x))\\ 2 Bt 



(x) 



\\f(y)\\ 2 y p(t,x,y)vol(dy). (22) 



Since for any h G L 1 (M), the function 



M 



h(y)p(t,z,y)vo\(dy) 



M 



is in C°°(M) (see Theorem 7.19 in [11]), we can use (14) with w = v 2 to 
deduce that for any compact K C M one has 



sup ( E 



i { *<c( 2) }e 2/ ^ 2(i?s(2 H / ||/(y)|i;p(*,^y)vol(dy) 



< oo, 



so that, in view of (22), the assignment 

^E[i {t<ax)} y t y/r 1 f(B t (x))] 



X H 



defines an element of F\_°o (M,E), and (21) is implied by the path integral 
formula from Theorem 2.9. ■ 

Next, we are going to deduce a finite propagation speed result, which will 
be used later on to prove that the compactly supported elements of D(H V ) 
are an operator core for Hy under geodesic completeness. The essential 
observation is that finite speed of propagation is always implied by a Davies- 
Gaffney type inequality, through a Paley- Wiener type theorem [3]. As we 
have already remarked in the introduction, we have borrowed this method 
from [13]. 
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Proposition 2.13. Let M be geodesically complete. 

(a) If V is bounded, then there is a constant D > such that for all open 
sets Ui, U 2 C M , all f ± , f 2 G r L 2(M, E) with supp(/j) C Uj and allt > one 
has 

\(e~ tHv h, f2)\ < e D V d ^^) 2 /^||M|||/ 2 ||. (23) 

(b) Let V be as in Theorem 2.6 and assume Hy > 0. Then for any compactly 
supported f G F L 2(M, E) and any t > 0, the section cos{t\/ H v ) f has a 
compact support. 

Proof, (a) Under the assumption that V is bounded and nonnegative, we 
are going to prove (23) with D = 0, which of course proves the assertion. 
To this end, we are going to use the well-known exponential-weight method, 
that goes back to [10] (see also [3]): Let q: M — > R be a bounded Lipschitz 
function with \dq\ < C a.e. in M. For any / G Tq^{M,E), Lemma 2.8 and 
the Sobolev product rule 



V( e 29 e-^/) = de 2q <g> e~ tHv f + e 2q Ve~ tHv f (24) 



imply 



A|| e «e-^/|| 2 

= -Re (^Ve~ tHv f^e- tHv f) - 2{Ve~ tHv f^e- tHv f) 
= -2Re (e q Ve- tHv f,e q dq®e- tHv f)-2\\e q Ve- tHv f\\ 2 
-2(Ve- tHv f,e 2q e~ tHv f). (25) 

Using Cauchy-Schwarz on the fibers for the first term and V > for the last 
term, the latter expression can be estimated by 

<2 / e^)||Ve-^/(x)|| x e^ \\dq(x)\\ x ||e-^/(x)|| x vol(dx) 

J M 

-2||e 9 Ve-^/|| 2 , (26) 
which, using XY < X 2 + Y 2 /4, is 

<i||e^|dg|e-^/|| 2 <^pe-^/|| 2 . (27) 

Thus, setting Sf, q {t) := ||e 9 e~'^ v /|| 2 , putting everything together and using 
Gronwall, we arrive at 

* f S) < e C2 * /2 ^(0). (28) 
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Now let U 1 ,U 2 be disjoint, let / G r C g°(M, E) with supp(/) C U 2 , and let 
a, b > 0. Then the function 

q : = ^min(6, ad(«,C/ 2 )) 
is bounded and Lipschitz with |dg| < a/2 a.e. in M and (28) implies 
\\lu^- tHv ff 

< e -min(6,ad([/ 1 ,C7 2 )) e 4* ( ^^ ) 

= e -min(b,ad([/ 1 ,C/ 2 )) e ^ /" || y ^ || 2 e min(b,ad(x,£/ 2 )) yol ( d;r ) 

•/t/ 2 

_ e -min(6,ad((7i,i7 2 )) e 4 1 \\f\\ 2 , (29) 



so that (by choosing a, 6 appropriately) 

HW^/H <e- d ^^) 2 /( 4t ) ll/H, (30) 
which carries over to f 2 by a density argument. Finally, we have 

|<e-^A,/ 2 >| = \{hAu^ tHv h)\ < e-^^^IIMHIMI (31) 

by Cauchy-Schwarz and (30), and everything is proved. 

(b) It is sufficient to prove that for any Uj, fj as in (a) and any < s < 
d(C/i, U 2 ) one has 

(cos(s v ^v)/ 1 ,/ 2 ) = 0. (32) 

Indeed, the latter implies that if supp(/) C K r (x) for some r > 0, x G M, 
then for any t > one has 

supp (cos (ty/W^j fj C K r+t (x), (33) 

and the latter set is compact by the geodesic completeness of M. It remains 
to prove (32). 

If V is bounded and Hy > 0, then (32) follows directly from (a): Indeed, 
one can use the same arguments as those in the proof of theorem 3.4 in [3] to 
see this. Essentially, one has to use a variant of the Paley- Wiener theorem, 
which has to be applied to an appropriately rescaled version of the analytic 
function z h-» (e~ zIiv /i, / 2 ), Rez > 0. 

Next, we assume that V is locally integrable and bounded from below 
with H v > 0. Then putting V n := min(V,n) for n G N (in the sense of the 
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fiberwise spectral calculus of E) we find by the above that (32) is satisfied 
for V replaced with V n , but monotone convergence of quadratic forms (see 
the proof of theorem 2.11 in [15]) gives Hy n — > Hy as n — > oo in the strong 
resolvent sense, which implies (32). 

Finally, if V is as in Theorem 2.14 and H v > 0, let us set V n := 
max(— n, V). Then each V n is locally integrable and bounded from below 
with Hy n > and again everything follows from the above and monotone 
convergence of quadratic forms (this is also included in the proof of theorem 
2.11 in [15]). ■ 

Now we are in the position to prove the main result of this paper: 

Theorem 2.14. Let M be geodesically complete, let \V\ G L^ oc (M) and as- 
sume that V has a decomposition V — V\ — V 2 into potentials Vj > with 
\V 2 \ G /C(M). Then r C g°(M, E) is an operator core for H v and one has 

D(H V )= {/|/,(V t V + y)/er L2 (M, J B)}. (34) 

Proof. We have to prove that r C go(M, E) is dense in D(H V ) with respect to 
the graph norm \\»\\ H ■ This will be proven in four steps: 
(I) Ifx e Cg°(M) and f G D(H V ), then X f e D(H V ) and 

H v (xf) = xHvf - V (dx)t / - l(A X )f. (35) 

Here, (d%)" denotes the vector field corresponding to the 1-form dx (with 
respect to the underyling Riemannian metric). 

Proof. We first note that the Sobolev product rule 

V(x/) = (d X ) ® / + xV/ (36) 

(which is applicable in view of (8) and (17)) shows that xf is m ^(qh v ), 
so that in order to prove xf G D(Hy), it is sufficient to construct a u G 
T L2 (M,E) such that 

QH v (xf,h) = {u,h} (37) 

for all h G D(q Hv ), where then H v (xf) is given by u. To this end, we 
calculate 

QH v (xf, h ) 

= \(y{xfivh) + (y{xflh) 

= ^(V/, V( X h)) - i(V/, (d X ) ®h) + ^((dx) ® /, V/i) + (Vf, X h) 

= (H v f, x h) - (v (dx) »M> + i((dtd x )/,/i>, 
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where we have used (36) in the second equality, and / G D(H V ) together 
with an integration by parts formula (Lemma 8.8 in [1]) and the Sobolev 
product rule 

® /) = (<!*«)/ -V a ,/ 

for (sufficiently) smooth 1-forms a in the third equality. In particular, we 
found a candidate u in (37) and it has the desired form as in (35). |B 

(II) The space 

D°(H V ) := D(H V ) fl j/ / has a compact support^ 
is dense in D(Hy) with respect to ||«||h v . 

Proof. By adding a constant, we can assume that H v > 0. But then the 
result readily follows from combining Proposition 2.13 with Lemma B.l. |l 

(III) The space 

D ooMx (H v ):=D (H v )nr l& (M 1 E) 
is dense in D°(Hy) with respect to ||»||h v - 

Proof. Let / G D°(H V ) and take r > 0, y G M with supp(/) C K r (y). 
Furthermore, pick a \ G Cg°(M) with x = 1 in K r+1 (y) and set f t := X e ~ tHy f 
for any t > 0. Then Proposition 2.11 implies f t G D^ oloc (iJv r ) and clearly 
Wft ~ /II — > as t — )> 0+. Furthermore, (I) implies H v (xf) — H v f and also 

^v(/t - /) = xHve~ tHv f - V (dx) »e-^/ - i(Ax)e"^/ - x^v/ 
+ V (dx)8 /+^(A X )/. 

Now it is easily seen that ||#y(/t - /)|| ->■ as t ->■ 0+. §■ 

(IV) r C g°(M, E 1 ) zs dense m D^ oloc (ifv p ) respect to \\*\\ H and one has 

d(h v )= {/|/,(v t v + y)/Gr L2 (M,£;)}. (38) 



Proof. Let / G loc (i?v). By Lemma 2.8 and the self-adjointness of H v we 
have H v C -£^> min , but it is well-known that (see for example p. 644 in [1]) 

D(^,J = {/| /, (V f V + V)f G T L2 (M, E)}. 

In particular D^ loc (// y ) C D(#* min ) implies iu := V f V/ + V7 G r L2 (M,£). 
As / is locally bounded with a compact support, one also has V f G T L 2 (M, i?), 
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so that V*V/ = w — Vf G T L 2(M, E). But now the assertion follows directly 
from Proposition A.l, which is in fact a local result (and which again heavily 
uses that / is locally bounded with a compact support). 

Finally, (38) simply follows from the essential self-adjointness of H Vmin , which 
follows from (II) and the by now proven fact that r C g°(M, E) is an operator 
core for Hy. ■ 

We immediately get: 

Corollary 2.15. Theorem 1.1 holds, that is, under the assumptions of The- 
orem 2.14) the operator V^V /2 + V is essentially self-adjoint on Tc^(M, E) , 
and its closure is semibounded from below. 

Proof. Combining Theorem 2.14 with Lemma 2.8 immeadiately gives Hy >m i n = 

H V,min = H V . ■ 



3 Application to Hydrogen type problems on 
Riemannian 3-manifolds 

In this section, we shall explain a typical application of Theorem 1.1: The es- 
sential self-adjointness of nonrelativistic Hamiltonians corresponding to Hy- 
drogen type atoms, with the electron's spin is taken into account. To this 
end, let us first explain what the analogues of the Coulomb potential and 
the Pauli operator are in a general curved setting. Here, we are going to 
follow [16] closely. 

Throughout Section 3, we will assume that M is a smooth con- 
nected Riemannian 3-manifold without boundary. 

Firstly, we want to point out that "nonparabolicity" is the appropriate setting 
that admits natural analogues of the Coulomb potential: 

Definition 3.1. The Riemannian manifold M is called nonparabolic, if one 
has 



pit, x, y)dt < oo for some (any) x,y G M with x ^ y. 



Then 



poo 

G: M x M — >(0,oo], G(x,y):= / p(t,x,y)dt 

Jo 

is called the Coulomb potential on M. 



16 



It should be noted that nonparabolicity always implies noncompactness. 
The essential point for the interpretation of G as the Coulomb potential is 
that M is nonparabolic, if and only if M admits a positive Green's function, 
and then G is the minimal positive Green's function (see [16] and the refer- 
ences therein for these facts). The following criterion can be easily deduced 
from Theorem 2.3: 

Lemma 3.2. Assume that there is a C > such that for allt > one has 

sup p(t,x,x) < Ct~ 3/2 . (39) 

xeM 

Then M is nonparabolic and there is a C > with 

G(x, y) < for all x,yeM. (40) 

d(x,y) 

Next, we will explain the natural analogues of the Pauli-operator in our 
general setting. To this end, we give ourselves a Pauli-Dirac structure (c, V) 
on M in the sense of [16], that is, with a smooth Hermitian vector bundle 
E -> M with rank£ = 2, 

c: T*M — >■ End(E) 
is a Clifford multiplication 4 and V is a Clifford connection 5 with respect to c. 

Remark 3.3. The existence of a Pauli-Dirac structure on M is a topological 
restriction, namely, M admits a Pauli-Dirac structure, if and only if M is a 
spin c manifold. This fact has also been explained in [16]. 

The Pauli-Dirac operator @>(c, V) with respect to (c, V) is defined by 

0(c,V) :=coV: Vc~(M,E) — > r c =«,(M, £), 

which is a linear first order differential operator with @(c, V)^ = 3>{c, V). If 
(ej) is some smooth local orthonormal frame for TM, then one has £^(c, V) = 
c ( e j)V ej - Furthermore, @(c, V) 2 is a generalized Laplacian on M which 
is given by the following Lichnerowicz formula: 

4 A Clifford multiplication c is a morphism of smooth vector bundles such that for all 
a e 1 (M) one has 

c(a) = —c(a)*, c(a)*c(a) = \a\ 2 . 

5 A Clifford connection is a Hermitian connection with the following property: for all 
a e O^M) and all X e r C oo (M, TM), V € r C oo(M,^) one has 

Vx(c(a)V) = c(V^ M a)^ + c(a)Vjc^. 
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Lemma 3.4. The differential form tr[V 2 ]/i G Vt 2 (M) is real-valued and 
closed, and one has 

9{c,V? = VtV + ±scai(.)l + Ij> [V 2 ] (e„ e,) C (e*) C (e*). (41) 

i<j 

The last lemma makes it plausible (see also Remark 3.6 (b) below) to call 
<^(c, V) := ^(c, V) 2 the Pauli-Dirac operator with respect to (c, V). 

Clearly, if one has (39), then G(», y) exists and is locally square integrable 
for any y G M, and for any such y and k, > one can consider the operator 

^(c,V;«, 2 /):=^(c,V)-«G(«, 2 /)l 

in I\2(M, £") with domain of definition Tc°°(M, i?), which gives rise to a 
symmetric operator. Let us furthermore define 

5(c,V):= / iscalWl + i^trfV^e^e,)^*)^) \ol(d:r) 
G [0,oo], 

where 1 1 1 * 1 1 1^ stands for the fiberwise Hilbert-Schmidt norm. With these 
preparations, Theorem 1.1 has the following important consequence: 

Theorem 3.5. Assume that M is geodesically complete with (39), and let 
S(c, V) < oo, k > 0, and y G M. Then iJ(c, V; n,y) is essentially self- 
adjoint and its closure H(c, V; n,y) is semibounded from below. 

Proof. Using ||«|| x < |||«||| x , (40) and Proposition 2.4, one easily checks that 
Theorem 1.1 can be applied with 

V := iscal(.)l + i^tr [V 2 ] (e 4 , ej )c(e*)c(e*) - «(?(•, 

i<j 

which proves the claim. ■ 

Remark 3.6. (a) A variant of Theorem 3.5 has been deduced in [16] with 
completely different methods, namely, using results of [1] (which rely on pure 
PDE methods). 

(b) In the situation of Theorem 3.5, the operator H(c,V; n,y) can be in- 
terpreted [16] as the nonrelativistic Hamiltonian corresponding to an atom 
with one electron and a nucleus with ~ k protons, where the electron's 
spin has been taken into account and the nucleus is considered as fixed in 
y with respect to the electron. Here, in view of Lemma 3.4, the underlying 
magnetic field is given by tr[V 2 ]/i G Q 2 (M). In particular, the assumption 
S(c, V) < oo is reasonable from the physics point of view, for it corresponds 
in a certain sense to a "finite magnetic self-energy" (it is essential for this 
interpretation to take the Hilbert-Schmidt norm in the definition of S(c, V)). 
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A Friedrichs mollifiers 

We record the following result on Friedrichs mollifiers here. Let M be 
a smooth connected Riemannian manifold without boundary, E — > M a 
smooth Hermitian vector bundle, V a Hermitian covariant derivative in E, 
and V : M ->■ End(E) a potential. 

Proposition A.l. Let \V\ G Lf oc (M) and assume that f G r L «. (M,E) is 
compactly supported with V^V/ G r L 2 (M, E) in the sense of distributions. 
Then there is a sequence (f n )nen C r c °°(M, E) such that 

lim ||/ B -/|| = 0, 

lim HV+V/n-VtV/H =0, 

lim \\Vf n -Vf\\=Q. 

n— s-oo 

Remark A. 2. Note that one indeed has / G T L 2(M,E), which follows from 
/ G i\oo (M, E) and the fact that / has a compact support. Furthermore, 
V f V/ Gr L2 (M, E) follows from V f V/ G r L 2 oc (M, E) and the fact that V f V/ 
has a compact support. 

Proof of Proposition A.l. Since most of the arguments should be well-known, 
we only sketch the proof. Let m := dim M and let d be the fiber dimension of 
E. Since / is compactly supported, we can use a partition of unity argument 
to assume that / is supported in a coordinate domain U C M (identified with 
an open subset of R m ) such that there is a smooth orthonormal frame for E 
over U, and we denote the components of / in this frame with f^\ . . . , f( d \ 
Now take some < j r G Cg°(M m ) with j(z) = for \z\ > 1 and 

f j(z)dz = l. 
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For r > let j r G Cg°(M m ) be given by j r (z) = r^jir^z). Let r > be 
small enough in the following such that the functions 

a— ► / j r (x-y)f {i) (y)dy, i = l,...d, (42) 

define an element 

f r er c ~(u,E)cr c ~(M,E). 

Since the sections f r — f and V^V/ r — V^V/ are compactly supported, the 
convergence 

lim ||/,- /||=0 (43) 

follows from Lemma 5.13 (ii) in [1], and 

lim ||V f V/ r - V f V/|| =0 

follows from the L^ oc - version of Proposition 5.14 in [1], which can be proven 
with analogous arguments. Note that so far we have only used that / is 
locally square integrable with a compact support. 

The local boundedness assumption on / comes into play as follows: Namely, 
this assumption combined with the compact support assumption implies that 
/ is actually bounded and so (42) implies 

ll/r(*)L<ll/IL forallx,r. (44) 

Since (in view of (43)) we may assume that f r —>f a.e. in M, and since f r 
has a compact support, the required convergence 

lim jVfr - Vf\\ = 
now follows from (44) and dominated convergence. ■ 



B Finite speed of propagation 

The following lemma is usually referred to as Chernoff 's finite speed of prop- 
agation method [2]. Let M be a smooth connected Riemannian manifold 
without boundary, and let E — > M be a smooth Hermitian vector bundle. 

Lemma B.l. Let S be a self-adjoint nonnegative operator in r L 2(M, E). 
Assume furthermore that D°(S), the compactly supported elements of'D(S), 
are dense in T L 2(M, E) and that for any f e D°(S') and any t > 0, the section 
cos(ty/~S)f has a compact support. Then D°(S') is an operator core for S. 

Proof. The proof is a straightforward generalisation of the proof of Theorem 3 
in [13]. ■ 
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C Path ordered exponentials 



In the following lemma, we collect some known facts about path ordered 
exponentials for the convenience of the reader: 

Lemma C.l. Let ffl be a finite dimensional Hilbert space, let T G (0, oo] 
and let F G L 1 1 oc ([0, T), Jz? Then the following assertions hold: 

(a) There is a unique weak (= AC locy ) solution Y : [0, T) — > Jzf(J^) of the 
ordinary initial value problem 

±Y(t) = Y(t)F(t), Y(0) = 1. (45) 

(b) For any < t < T one has 

oo „ 

Y(t) = 1 + / F ( s • • • F ( s k)d Sl . . . ds k . (46) 

k=l ^0< S1 <---<s k <t 



(c) If F(») is Hermitian a.e. in [0, T) and if there exists a real-valued func- 
tion c G Lj L oc [0, T) such that for all v G it holds that 

(F(»)v, v)jg, < c(«) \\v\\^, a.e. in[0,T), 

then one has 

\\Y{t)\\^ < eti< s)ds for allO <t<T. 
Proof. See [7] and the Appendix C of [15]. ■ 
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